Dynamics of cubic-tetragonal phase transition in KNDO3 perovskite 
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The low-energy part of the vibration spectrum in KNbC>3 was studied by cold neutron inelastic 
scattering in the cubic phase. In addition to acoustic phonons, we observe strong diffuse scattering, 
which consists of two components. The first one is quasi-static and has a temperature-independent 
intensity. The second component appears as quasi-elastic scattering in the neutron spectrum in- 
dicating a dynamic origin. From analysis of the inelastic data we conclude that the quasi-elastic 
component and the acoustic phonon are mutually coupled. The susceptibility associated with the 
quasi-elastic component grows as the temperature approaches Tc- 

PACS numbers: 77.80.-e, 61.12.-q, 63.50.+X, 64.60.-i 



ABO3 perovskites form a class of important materi- 
als, in part because of potential technical applications 
but also as fundamental interest in the physics of phase 
transitions [1, 2]. At sufficiently high temperatures many 
of these perovskites have 0\ cubic symmetry and struc- 
tural phase transitions can take place as the tempera- 
ture is lowered. Well-known examples are e.g. the cubic- 
tetragonal phase transition in SrTi03 (Tc w 105 K) or in 
BaTiO-3 (T c « 425 K) (for a review see Ref. [3]). There 
are, however, ABO3 perovskites which were less stud- 
ied. An example is the first-order cubic-tetragonal phase 
transition in KNb03 which occurs at Tc ss 683 K when 
cooling the crystal from above the transition tempera- 
ture [4]. 

The mechanism of the cubic-tetragonal (CT) phase 
transition in KNb03 is still controversial. Whereas well- 
defined soft phonon modes with frequency varying with 
temperature have been detected in many materials close 
to Tc [1 - 3], only an over-damped excitation has been 
observed in cubic KNb03 with neutron scattering and 
it was suggested that the nature of the C-T phase tran- 
sition in that compound is similar to the displacive C-T 
transition in BaTi03 [5, 6]. On the other hand, two coex- 
isting and essentially uncoupled modes are inferred from 
analysis of optical data in the cubic phase of KNb03 : a 
relaxation mode and a soft phonon, with the relaxation 
process driving the C-T phase transition [7]. 

We re-investigated the low energy part of the vibration 
spectrum in KNb03 under improved resolution condi- 



tions first to try to elucidate the mechanism of the phase 
transition in this crystal and second to check whether 
the diffuse scattering found in Ref. [8] is of static or dy- 
namic origin. The inelastic cold-neutron scattering mea- 
surements reported here were performed with the three- 
axis spectrometer TASP, located at the neutron spal- 
lation source SINQ (Paul Scherrer Institute, Switzer- 
land). A large single crystal of KNb03 (~ 20 cm 3 , 
mosaic ~ 80') was mounted into an ILL-type furnace. 
To decrease the level of incoherent background the sam- 
ple holder was made from pure niobium. The crys- 
tal was aligned in the (h k 0) scattering plane. The 
measurements were performed in the temperature range 
727 K - 1030 K. The (002) reflection of pyrolytic graphite 
(PG) was used to monochromate and analyze the in- 
cident and scattered neutron beams, respectively. The 
spectrometer was operated in the constant final-energy 
mode with kj = 1.97A . A PG filter was used to re- 
move higher-order wavelengths. The horizontal collima- 
tion was 10'/ A - 80' - 80' - 80'. With that configu- 
ration the energy resolution at zero energy transfer is 
~ 0.4 meV. By monitoring the position and intensity of 
the (1, 1, 0) Bragg peak, the temperature of the cubic- 
tetragonal phase transition upon cooling was found at 
Tq = 684 ± 2 K, in close agreement with published data 
[4, 6]. The temperature of the sample was controlled by 
two thermocouples. The temperature gradient through 
the sample did not exceed 15 K. 

Before analyzing the inelastic neutron spectra quan- 
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two Lorcntzian components (see inset of Fig. 1). 
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FIG. 1: Profiles of elastic scans in the (OryO) direction at 
T = 760 K and T = 1030 K . Raw data are shown by circles. 
The dotted line is the result of fits as described in the text. 
The bold line shows the Lorentzian profile, the dashed-dotted 
line stands for the intense and narrow Bragg peak. The in- 
tensity is given in a logarithmic scale. The insert shows the 
difference of the elastic scans I(T = 760K)-I(T = 1030K) fit- 
ted with Eq.0 1 rlu corresponds to 1.57 A -1 . 



titatively, it is convenient to address the q-dependence 
of the elastic neutron response. Figure ^ shows repre- 
sentative elastic scans along the (l,l±q,0) direction at 
T=760 K and 1030 K, respectively. This intense and 
broad scattering is similar to the diffuse scattering ob- 
served in KNbC>3 by Guinier et al. [8] using X-rays and 
reflects the presence of atomic disorder in the perovskite 
cell. In KNbC>3 atomic disorder yields diffuse scattering 
along the [100] direction both in the X-ray and neutron 
diffraction patterns. Here we approximate the line-shape 
of the neutron diffuse scattering intensity by a Lorcntzian 
profile: 
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where q is the position in reciprocal space; k the inverse 
of the correlation length £ and Iq yields the integrated 
intensity. From a fit to the elastic data at T = 1030 K 
we obtain £ = 64 ± 6 A. It turns out that the shape and 
intensity of the diffuse scattering measured in the (2,0,0) 
Brillouin zone (BZ) does not depend on temperature (see 
insert of Fig. |2J . This is in agreement with the results 
of Rcf. [8] where the intensity of the diffuse scattering is 
found to be temperature independent in the cubic phase 
and to decrease abruptly by ~ 30% immediately below 
Tcr. On the other hand, it turns out that the inten- 
sity of the diffuse scattering measured along (l,l±g,0) 
slowly decreases when increasing the temperature from 
To- This suggests that in that BZ and for temperatures 
relatively close to Tq, the diffuse scattering consists of 
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FIG. 2: Neutron scattering spectrum from KNbC>3 at 1030 
K. Raw data are shown by open circles. The solid line is the 
result of fit as described in the text. The intensity is given 
in a logarithmic scale. The apparent width of the phonon 
peak is due to the resolution effects. The insert shows the 
temperature and the q dependences of the intensity of the 
central peak. 1 rlu=1.567 A -1 . 

We turn now to the analysis of the inelastic neutron 
scattering spectra. Figure shows an example of a 
constant-q scan taken at Q =(2, -0.1, 0) and T = 1030 K. 
The spectrum contains an inelastic peak at fiu)=4 meV 
from the transverse acoustic (TA) phonon and a narrow 
peak centered around zero energy transfer. To analyze 
the data quantitatively, we hence modeled the neutron 
scattering intensity I(Q,uj) in the following way: 



J(Q,w) = S(Q,w)®R(Q,w)+B 



(2) 



The symbol (8> stands for the 4-D convolution with the 
spectrometer resolution function R(Q,lu) [9]; B denotes 
the background level; S(Q, ui) is the neutron scattering 
function which is related to the the imaginary part of the 
dynamical susceptibility x"(Q,w) through 



S{Q,u) 



(3) 



with the temperature factor [n(u>) + 1] = [1 — 
exp(— w/T)] -1 . We approximate the central peak by a 
(5-function in energy 



S CP = A(q)S(u) 



(4) 



with A(q) given by Eq.^ The line-shape of the acoustic 
phonon is given by the usual damped-harmonic oscillator 
(DHO) 



xdho(q, w) = (fig - «7?^ - 



(5) 
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In Eq. [5] 7 9 is the damping and Q, q = 

Lu q = c • q [10] is the rcnormalizcd frequency of the acous- 
tic phonons. For small values of momentum transfers q, a 
linear dispersion for the acoustic phonon branch is a rea- 
sonable approximation and the phonon damping approx- 
imatively follows a c?q 2 -dependence [11]. The scattering 
function used to fit the neutron data then reads 



S(Q,lo) = Scp(Q,w) 



K^) + jj t 2 n / \ 



(6) 



where Q = q + t is the neutron scattering vector and r 
a reciprocal lattice vector; /i is the structure factor of 
the acoustic phonon. As shown in Fig[3 Eqs. 1-6 param- 
eterize the experimental data in the (2, 0, 0) BZ well. 
The central peak is resolution-limited and temperature- 
independent. The acoustic phonon branch has a stiffness 
c = 28 ± 1.3 meV-A 2 and the damping is small at low q, 
7 9 = dq 2 with d = 55 ± 6 meVA 2 (0.05< q < 0.15 (rlu)). 
In the temperature range 750 K< T <1030 K no qualita- 
tive change in the dispersion of the acoustic phonon was 
observed for data taken along (2, q, 0) (\q\ <0.15 (rlu)). 

On the contrary, for temperatures below T=1030 K an 
additional component is observed in the inelastic spectra 
for constant-q scans in the (1,1,0) BZ. For means of com- 
parison, Fig.[3]shows two representative neutron scatter- 
ing spectra measured in the (1, 1, 0) BZ at T=1030 K 
and T=727 K, respectively. At T = 1030 K the spec- 
trum consists of two components - a central peak (CP) 
and a phonon response around hu =4 meV. However, 
as the temperature is lowered to T=727 K, additional 
quasielastic scattering (QE) appears along (1, l±q, 0) 
direction. The intensity of this quasi-elastic scattering 
grows when approaching Tc- From the above discus- 
sion we conclude that in the (1, 1, 0) BZ, the inelastic 
neutron spectra consist of three contributions: a central 
peak, quasi-clastic and phonon scattering. To describe 
the quasi-elastic scattering we introduce a Debye-likc re- 
laxation function 



Xq-el(<l,u) 



X(0,T) 
1 + q 2 jn 2 



(i-iu/r q y 



(7) 



x(0, T) is the temperature dependent static susceptibil- 
ity; k the inverse of the correlation length and T q = 
To + Dq 2 . Taking into account the quasi-elastic scat- 
tering modifies the neutron cross-section to 



S C p(Q,uj) 
[n(u) + 1] 



(8) 
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However, the scattering function given in Eq. [S] fails in 
reproducing the experimental data in the (1, 1, 0) BZ 
for T<1030 K. For example Fig. [3] shows an inelastic 
spectrum at T=727 K and Q=(l, 1.075, 0) where a 
qualitative change in the phonon line-shape accompa- 
nied by a shift in the position of the phonon peak is ob- 
served. These two effects suggest that coupling between 
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FIG. 3: Neutron scattering spectrum from KNb0 3 at T=1030 
and 727 K, respectively. The solid and dashed line are the re- 
sults of fit as described in the text. To emphasize the QE 
component the intensity is given in a logarithmic scale. Note 
the pronounced change in the phonon line-shape at lower tem- 
perature. 



the quasi-elastic component and the acoustic phonons 
becomes important as the temperature approaches Tc. 
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FIG. 4: Observed and fitted inelastic neutron intensities taken 
in the (1, 1, 0) BZ at T=727 K. Fitted curves were obtained 
with Eq. I1UI Intensity is given in a logarithmic scale. 

The dynamical susceptibility for two coupled excita- 
tions was considered in details in Rcfs. [3, 12, 13] and is 
given by 

XCm{Q,U) = : T2 9 

1 - ^ X1X2 

where Xi = XiCQj^)) i = 1,2 are the dynamical suscep- 
tibilities of the uncoupled phonon and QE component, 
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respectively. In the following we take /j as real constants 
since in KNDO3 all the atoms are situated on centers 
of symmetry. The interaction term is A = A(q,u>) = 
(g r + iojgi) q 2 . Finally, the scattering function reads: 



S(Q,u) = S C p(Q,u J ) + 



[n(w) 



1] „ 

XCM ■ 



(10) 



In order to obtain a good agreement between Eg. 1101 and 
the neutron spectra, it was necessary to fit the com- 
plete set of data (0<q<0.2) taken at a given tempera- 
ture simultaneously. Figures 0] and [3] show the results 
of such calculations for T=727 K and T=760 K. We ob- 
tain g r = 20 ± 3 meV 2 • A 2 and g 4 = 95 ± 6 meV-A 2 . 
Introduction of a coupling between the QE and acoustic 
modes has two consequences. First, Eq. yields a bet- 
ter description of the line-shape of the inelastic neutron 
spectra. Second, Xcm(Q> w ) * s enhanced at low energy 
transfers. Further, we obtained To = 0.19 ± 0.05 meV 
and D = 44 ± 4 meV-A 2 for the damping of QE com- 
ponent. As discussed above, both the CP and the 
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FIG. 6: Temperature dependence of the susceptibility of the 
QE component. The solid line is a fit to the data as explained 
in the text. 



10000 

aj. 1000 
en 



A. KNb0 3 , T=760 K 


/ r *\\ 


• (1,1.05, 0) 






p (1,1.075,0) 




/ /rrVSv 


- (1,1.1, 0) 






A (1,1.15, 0) 




/ rn 




f jr Ik iV\n ^ 




ysA on. 










j!" it 4 — ■ ^xf^SgSBg 






2 2 


4 


6 



Energy (meV) 

FIG. 5: Observed and fitted inelastic neutron intensities taken 
in the (1, 1, 0) BZ at T=760 K. Fitted curves were obtained 
with Eq. 1101 Intensity is given in a logarithmic scale. 

line-shape of the acoustic phonons are temperature inde- 
pendent in the (2, 0, 0) BZ. Hence, to fit the data mea- 
sured in the (1, 1, 0) BZ as a function of temperature, we 
fixed the parameters of the CP and the acoustic phonons. 
The only parameter left to describe the temperature de- 
pendence of the neutron spectra is the susceptibility of 
QE scattering x(0,T). As shown in Fig. the intensity 
of the quasi-elastic component is a maximum close to 
Tc and decreases continuously with increasing temper- 
ature. The temperature dependence of x(0, T) follows 
approximately the Curie- Weiss law oc 1/(T — Tq) with 
To = 615 ± 23 K in good agreement with the value de- 
duced from dielectric measurements To = 633 ± 5.9 K 
(Ref. [14]), and T = 615 K (Ref. [15]). This suggests 
that the cubic-tetragonal transition in KNb03 is driven 



by the quasi-elastic relaxational excitation. The intensity 
of the quasi-elastic component is strong in the (1,1,0) 
zone and has a small intensity in the (2,0,0) Brillouin 
zone, which indicates that the relaxation mode is due 
to correlated atomic motion of optical character. How- 
ever, at all q and temperatures we did not observe that 
the relaxation mode evolves into an under-damped optic 
phonon branch. Thus, we conclude that QE scattering in 
KNb03 is not due to an usual overdamped soft phonon 
but is related to disorder in the lattice. 

To summarize, we measured the low-energy part of the 
vibration spectrum of KNb03 in the cubic phase with 
inelastic neutron scattering. We find a coexistence of a 
static and a quasi-elastic component. The static compo- 
nent appears to correspond with static disorder in the 
cubic cell and is temperature independent in agreement 
with X-rays results [5]. The quasi-elastic component is 
coupled with the acoustic phonon branch and its inten- 
sity follows the Curie- Weiss law well. 
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